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Abstract 

The basic features of the dynamics of open quantum systems, such as the dissipation of energy, 
the decay of coherences, the relaxation to an equilibrium or non-equilibrium stationary state, and 
the transport of excitations in complex structures are of central importance in many applications of 
quantum mechanics. The theoretical description, analysis and control of non-Markovian quantum 
processes play an important role in this context. While in a Markovian process an open system 
irretrievably loses information to its surroundings, non-Markovian processes feature a flow of in- 
formation from the environment back to the open system, which implies the presence of memory 
effects and represents the key property of non-Markovian quantum behavior. Here, we review 
recent ideas developing a general mathematical definition for non-Markoviantiy in the quantum 
regime and a measure for the degree of memory effects in the dynamics of open systems which 
are based on the exchange of information between system and environment. We further study the 
dynamical effects induced by the presence of system-environment correlations in the total initial 
state and design suitable methods to detect such correlations through local measurements on the 
open system. 
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I. INTRODUCTION 



The standard approach to the theoretical description and analysis of dissipation and deco- 
herence processes in open quantum systems presupposes a time evolution without memory. 
Employing the concept of a quantum Markov process which is given by a semigroup of 
completely positive dynamical maps one obtains a quantum Markovian master equation 
describing the time evolution of the reduced open system states with a generator in Lind- 
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blad form 121 • Within a microscopic approach quantum Markovian master equations 
are usually obtained by means of the Born-Markov approximation which assumes a weak 
system-environment coupling and several further, mostly rather drastic approximations. 
However, in many processes occurring in nature these approximations are not applicable, a 
situation which occurs, in particular, in the cases of strong system-environment couplings, 
structured and finite reservoirs, low temperatures, as well as in the presence of large initial 
system-environment correlations jsj]. In the case of substantial quantitative and qualitative 
deviations from the dynamics of a quantum Markov process one often speaks of a non- 
Markovian process, implying that the dynamics is governed by significant memory effects. 
However, a consistent general theory of quantum non-Markovianity does not exist and even 
the very definition of non-Markovianity is a highly topical issue. Very recently important 
steps towards a general theory of non-Markoviantiy have been made which try to rigorously 
define the border between Markovian and non-Markovian quantum dynamics and to de- 
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velop quantitative measures for the degree of memory effects in open systems [4-6]. These 
approaches provide general mathematical characterizations of quantum non-Markovianity 
which are independent from any specific representation or approximation of the dynamics, 
e.g. in terms of a perturbative master equation. 

The key questions to be studied here are, how can one mathematically define and quantify 
non-Markovian behavior in the quantum regime, how do quantum memory effects manifest 
themselves in the dynamical behavior of complex open quantum systems, and how can 
such effects be uniquely identified? The answer to these questions is of great relevance for 
the design of appropriate schemes allowing the experimental detection and measurement 
of non-Markovianity. Here, we review some recent ideas and concepts which characterize 
non-Markovian quantum behavior by means of t 
an open quantum system and its environment 



re information which is exchanged between 
The gain or loss of information can 



be quantified through the dynamics of the trace distance between a pair of quantum states 
of the open system. It is known that this distance measure for quantum states can be 



interpreted as a measure for the distinguishability of the states [8HlO|. Markovian processes 
tend to continuously reduce the trace distance and, hence, the distinguishability between 
any pair of physical states, which means that there is a flow of information from the open 
system to its environment. In view of this interpretation the characteristic feature of a non- 
Markovian quantum process is the increase of the distinguishability, i.e., a reversed flow of 
information from the environment back to the open system. Memory effects thus emerge 
through a recycling of information such that earlier states of the open system influence its 
future states. 

We first recapitulate in Sec. [Til all necessary concepts from the quantum theory of open 
systems, such as completely positive quantum dynamical maps, dynamical semigroups, Lind- 
blad generators, quantum master equations and the notion of the divisibility of dynamical 
maps. The precise definition of non-Markovian quantum dynamics and the corresponding 
measure for the degree of memory effects will then be discussed in detail in Sec. IHII This 
section also contains a summary of the most important mathematical and physical features 
of the trace distance. In Sec. |IV]we study a further important aspect, namely the dynami- 
cal influence of correlations in the initial system-environment state. We will derive general 
inequalities expressing this influence and develop experimentally realizable schemes to lo- 
cally detect such correlations. Finally, some conclusions and further developments will be 
discussed in Sec. [V] 



II. OPEN QUANTUM SYSTEMS: BASIC NOTIONS AND CONCEPTS 

We introduce and discuss some of the most important notions and general concepts of 
the quantum theory of open systems. The presentation is of course not an exhaustive 
review, but we concentrate on those topics that are needed for the discussions of quantum 
non-Markovianity in the following sections. 
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A. Microscopic representation of open systems 



An open quantum system S is a quantum system which is coupled to another quantum 
system E, its environment. Thus, S can be regarded as subsystem of the total system S + E 
consisting of open system plus environment. Denoting the Hilbert spaces of S and E by %$ 
and He, respectively, the Hilbert space of the total system is given by the tensor product 
space 

u = u s ®u E . (1) 

The physical states of the total system are described by density matrices p, representing 
positive trace class operators on H with unit trace. This means that p > 0, which implies 
that p is self-adjoint with nonnegative eigenvalues, and trp = 1. The corresponding reduced 
states of subsystems S and E are given by partial traces over He and Hs, respectively, 

ps = tr E p, Pe = tr s p. (2) 

In the following we denote the convex set of physical states pertaining to some Hilbert space 
H by S(H). 

We will assume here that the total system S + E is closed and follows a unitary dynamics 
described by some unitary time evolution operator 

U(t)=exp[-iHt] (h=l) (3) 

with a Hamiltonian of the most general form: 

H = H s ®I E + Is®H E + H z , (4) 

where H$ and He are the self-Hamiltonians of system and environment, respectively, and 
Hi is any interaction Hamiltonian. Thus, the time dependence of the total system states is 
given by the von Neumann equation, 

j t p(t) = -i[H,p(t)}, (5) 

with the formal solution 

p(t) = U(t)p{V)U\t). (6) 

It turns out that in most cases of practical relevance a full solution of the equations of 
motion on the microscopic level is impossible. Thus, one of the central goals of the quantum 



theory of open systems [3| is the development of an effective and analytically or numerically 
feasible formulation of the dynamical behavior of a suitably defined reduced set of variables 
forming the open subsystem S. Given a certain split of the total system into open system S 
and environment E, on tries to derive effective equations of motion for the time dependence 
of the reduced system state ps{t) through an elimination of the environmental variables from 
the dynamical equations. The main aim is to develop efficient descriptions for a wide class 
of physical problems and phenomena, such as the dissipation and damping of populations, 
the relaxation to a thermal equilibrium state, the emergence of non-equilibrium stationary 
states, the suppression or destruction of quantum coherences, the description of the quantum 
transport of excitations in complex systems, and the dynamics of quantum correlations and 
entanglement. 



B. Quantum dynamical maps 

Quantum dynamical maps represent a key concept in the theory of open systems. To 
introduce this concept we presuppose (i) that the dynamics of the total system is given by 
a unitary time evolution ([6]), and (ii) that system S and environment E are statistically 
independent at the initial time, i.e., that the initial total system state represents a tensor 
product state 

p(0)= Ps (0)®p E (0). (7) 

While condition (i) can be relaxed to include more general cases, condition (ii) is crucial 
for the following considerations. We will discuss in Sec. IIVI the case of initially correlated 
system-environment states. On the basis of these assumptions the open system state at time 
t > can be written as follows, 

p s (t) = tr E {U(t) Ps (0)®p E (0)rf(t)} . (8) 

In the following we consider the initial environmental state ps(0) to be fixed. For each fixed 
t > Eq. (jSJ) is then seen to represent a linear map 

$(t,0): S(Hs)^S(H s ) (9) 

on the open system's state space S(T-Ls) which maps any initial open system state to the 
corresponding open system state at time t: 

Ps(P)^Ps(t)=$(t,0)p s (p). (10) 
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This is the quantum dynamical map corresponding to time t. It is easy to check that it 
preserves the Hermiticity and the trace of operators, i.e., we have 

tisMt,0)A}=ti s {A}, (11) 

and 

[$(t,0)A] f = <5>(t,0)Al (12) 
Moreover, 0) is a positive map, i.e., it maps positive operators to positive operators, 

A>0 =>- $(£, 0)A > 0. (13) 

An important further property of dynamical maps is that they are not only positive but also 
completely positive. Such maps are also known as trace preserving quantum operations or 



quantum channels in quantum information and communication theory [ll|. We recall that 
a linear map $ is completely positive if and only if it admits a Kraus representation 12), 
i.e., if there are operators fli on the underlying Hilbert space Us such that 

®A = J2 n i m l ( 14 ) 

i 

Such a map is trace preserving if and only if the normalization condition 

£ntn = i 5 (is) 

i 

holds. The original definition of complete positivity, which is equivalent to the existence of 
a Kraus representation, is the following. Consider for any n = 1,2, . . . the tensor product 
Hs ® C n , describing the Hilbert space of S combined with an n- level system, and the 
corresponding linear tensor extension of $ defined by ($ (g I n ){A <8> B) = (&A) Cg B. The 
map $ (g) I n thus describes an operation which acts only on the first factor of the composite 
system, and leaves unchanged the second factor. The map $ is then defined to be n-positive 
if $ Cg) I n is a positive map, and completely positive if $ <g> /„ is a positive map for all n. 
We note that positivity is equivalent to 1-positivity, and that for a Hilbert space with finite 



13|- 



dimension Ng = dimT-Ls complete positivity is equivalent to ^-positivity 

If we now allow the time parameter t to vary (keeping fixed the initial environmental 
state pe{0)), we get a one-parameter family of dynamical maps, 

{$(t,0) |t>0,$(0,0) = /}, (16) 



which contains the complete information on the dynamical evolution of all possible initial 
system states. Thus, formally speaking a quantum process of an open system is given by 
such a one-parameter family of completely positive and trace preserving (CPT) quantum 
dynamical maps. 

As a simple example we consider the decay of a two-state system into a bosonic reservoir 
. The total Hamiltonian of the model is given by Eq. (J4]) with the system Hamiltonian 

H s = u a + a-, (17) 

describing a two-state system (qubit) with ground state |0), excited state |1) and transition 
frequency uq, where a + = |1)(0| and cr_ = |0)(1| are the raising and lowering operators of 
the qubit. The Hamiltonian of the environment, 

H E = J2u k b{b k , (18) 

k 

represents a reservoir of harmonic oscillators with creation and annihilation operators b\ and 
bk satisfying Bosonic commutation relations b k ,\ = 5kk'- The interaction Hamiltonian 
takes the form 

Hi = ^{g k a + ®b k + gla_®b ] ^ . (19) 

k 

The model thus describes for example the coupling of the qubit to a reservoir of electromag- 
netic field modes labelled by the index k, with corresponding frequencies and coupling 
constants g^. Since we are using the rotating wave approximation in the interaction Hamil- 
tonian, the total number of excitations in the system, 

N = a+a-+J2 b kh, (20) 
k 

is a conserved quantity. The model therefore allows to derive an analytical expression for 
the dynamical map ([TO]) . Assuming the environment to be in the vacuum state |0) initially 
one finds: 

Pll (t) = |G(t)| 2 pn(0), (21) 
Poo(t) = Poo(0) + (1 - \G(t)\ 2 )p n (0), (22) 
p 10 (t) = G(t)p w (0), (23) 
Poi(t) = G*(i)p i(0), (24) 



where the pijit) = (i\ps(t)\j) denote the matrix elements of ps(t). The function G(t) 
introduced here is defined to be the solution of the integro-differential equation 



^G(t) = -J*dt 1 f(t-ti)G(t 1 ) (25) 



corresponding to the initial condition G(0) = 1, where the kernel fit — t\) represents a 
certain two-point correlation function, 

f{t-t x ) = {0|5(t)^ t (t 1 )|0>e^ o{t " tl) 

= J2\9k\ 2 e i{ul0 ~ ulk){t ~ tl \ (26) 
k 

of the environmental operators 

B(t) = Y,9kb h e-*"' t . (27) 
k 

These results hold for a generic environmental spectral density and the corresponding two- 
point correlation function. Taking, for example, a Lorentzian spectral density in resonance 
with the transition frequency of the qubit we find an exponential two-point correlation 
function 

f(r) = ^7oAe- AM , (28) 

where 70 describes the strength of the system-environment coupling and A the spectral width 
which is related to the environmental correlation time by te = A . Solving Eq. ( 125]) with 
this correlation function we find 



G(t) = e~ xt / 2 

where d = a/A 2 — 270A. 



, , dt\ A . . ( dt 
cosh(- +-smh - 



(29) 



C. Completely positive semigroups 

The simplest example of a quantum process is provided by a semigroup of completely 
positive dynamical maps, which is often considered as prototypical example of a quantum 
Markov process. In this case one assumes that the set ffl6|) has the additional property 



$(t,0)$(s,0) = $(t + s,0) 
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(30) 



for all t, s > and, hence, has the structure of a semigroup. Under very general mathematical 
conditions such a semigroup has an infinitesimal generator C which allows us to write 

$(t,0)=exp[£t]. (31) 

Accordingly, the reduced system state ps(t) obeys the master equation 

j t p s (t) = Cp s {t). (32) 

The complete positivity of the semigroup leads to important statements on the general 
structure of the generator. The famous Gorini-Kossakowski-Sudarshan-Lindblad theorem 



[i 



2j states that C is the generator of a semigroup of completely positive quantum dynamical 
maps if and only if it has the following form, 



C Ps = -i [H s , p s ] + A ^A - 2 {A*** Ps} 



(33) 



where Hs is a system Hamiltonian (which need not coincide with the system Hamiltonian 
Hs in the microscopic Hamiltonian (jlj)), the Ai are arbitrary system operators, often called 
Lindblad operators, describing the various decay modes of the system, and the 7$ are corre- 
sponding decay rates. This theorem has many far-reaching consequences and is extremely 
useful, in particular in phenomenological approaches since it guarantees a time evolution 
which is compatible with general physical principles for any master equation of the above 
structure. On the other hand, it is in general difficult to justify rigorously the assump- 
tion ( l30i) and to derive a quantum master equation of the form ( 132]) starting from a given 
system-environment model with a microscopic Hamiltonian (jlj). Such a derivation requires 
the validity of several approximations, the most important one being the so-called Markov 
approximation. This approximation presupposes a rapid decay of the two-point correlation 
functions of those environmental operators that describe the system-environment coupling 
H\. More precisely, if te describes the temporal width of these correlations and tr the 
relaxation or decoherence time of the system, the Markov approximation demands that 

t e <^t r . (34) 

This means that the environmental correlation time te is small compared to the open sys- 
tem's relaxation or decoherence time tr, i.e., that we have a separation of time scales, the 
environmental variables being the fast and the system variables being the slow variables. 



The Markov approximation is justified in many cases of physical interest. Typical examples 

of application are the weak coupling master equation, the quantum optical master equa- 

fl 

tion describing the interaction of radiation with matter [3], and the master equation for 
a test particle in a quantum gas [15]. However, large couplings or interactions with low- 
temperature reservoirs can lead to strong correlations resulting in long memory times and in 
a failure of the Markov approximation. Moreover, the standard Markov condition (1341) alone 
does not guarantee, in general, that the Markovian master equation provides a reasonable 
description of the dynamics, a situation which can occur, for example, for finite and/or 
structured reservoirs Q, Q ■ 



D. Time-local master equations 



There exists a whole bunch of different theoretical and numerical methods for the treat- 
ment of open quantum systems, beyond the assumption of a dynamical semigroup, such as 
projection operator techniques 



18 



19], influence functional and path integral techniques 

2- Here, 



20J, quantum Monte Carlo methods and stochastic wave function techniques 
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we concentrate on a specific approach which is particularly suited for our purpose and which 
describes the open system dynamics in terms of a time-local master equation. 

It is usually expected that the mathematical formulation of quantum processes describing 
effects of finite memory times in the system must necessarily involve equations of motion 
which are non-local in time. In fact, such a description is suggested by the Nakajima-Zwanzig 
projection operator technique which leads to an integro-differential equation for the reduced 



density matrix 18 



19] . However, even the presence of strong memory effects does not 



exclude the description of the dynamics in terms of a quantum master equation which is 
local in time, as may be seen from the following simple argument. According to Eq. (flOl) we 
have ps(t) = 0)ps(0). Assuming a smooth time-dependence we may differentiate this 
relation to get 

ri 

(35) 



-p s (t) = $(t,0)p s (0), 

where the dot indicates the time derivative of <£>(£, 0). To obtain a local master equation we 
invert the relation (TTU|) . expressing ps(0) in terms of ps(t), which yields 

d 



dt 



p s (t) = <i>(t,0)<5>- 1 (t,0)ps(t). 



(36) 
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Thus we see that the linear map JC(t) = $(t, 0)$ _1 (t, 0) represents a time-dependent gen- 
erator of the dynamics and we obtain a quantum master equation which is indeed local in 
time, providing a linear first-order differential equation for the open system state: 

j t p s (t) = K(t)p 8 (t). (37) 

We note that the above argument presupposes that the inverse $ _1 (t, 0) of the map $(t, 0) 
exists. It is possible that the inverse of $(t,0) and, hence, also the time-local generator 
K,{i) do not exist. Such a situation can indeed occur for very strong system-environment 
couplings (see below). However, for an analytic time- dependence the inverse of 0) and 



the generator JC(t) do exist apart from isolated singularities of the time axis [22]. It should 
also be emphasized that 0) denotes the inverse of $(t, 0) regarded as linear map acting 

on the space of operators of the reduced system. The important point is that this does not 
imply that $ _1 (t, 0) is required to be completely positive. In general, the inverse map is 
not only not completely positive, but even not positive. There exists a powerful method for 
the microscopic derivation of time-local master equations of the form ( 1371) which is known 
as time-convolutionless projection operator technique js], [23L 27 ] . This technique yields a 



systematic expansion of the generator of the master equation in terms of ordered cumulants 
and many examples have been treated with this method (28^32 ] . 

The generator JC(t) of the time-local master equation must of course preserve the Her- 
miticity and the trace, i.e., we have 

[IC^A^ = K{t)A\ (38) 

and 

tr 5 {!C(t)A} = 0. (39) 

From these requirements it follows that the generator must be of the following most general 
form, 



K(t)ps = -i[H s (t),p s ] 



+ !>(*) Ui(t)p s At(t) - \ {At(t)Mt),p s } 



(40) 



We see that the structure of the generator provides a natural generalization of the Lindblad 
structure, in which the Hamiltonian Hs(t), the Lindblad operators Ai(t) as well as the 
various decay rates ji(t) may dependent on time. 

11 



We briefly sketch the proof of (140]) for a finite dimensional open system Hilbert space, 
dim Hs = Ng. To this end, we consider a fixed time t and a fixed complete system of 
operators Fi on U s , i = 1,2,..., N 2 S , which are orthonormal with respect to the Hilbert- 
Schmidt scalar product, 

txsiFjFj} = 5 ir (41) 
Without loss of generality, we may further choose 



F "i = 7wl Is - (42) 



According to Lemma 2.3 of Ref. [1] any linear map K(t) satisfying Eqs. ( 1381) and ( 1391) can 
then be written as 



K(t) Ps = -i [H s (t), Ps ] 



Nl-l 



1 

+ £ % (t) f^f;--!^,^} 



(43) 



with a self-adjoint operator Hs(t) and a Hermitian matrix c(i) = (cjj(t)). Diagonalizing this 
matrix with the help of a unitary matrix u{t) = (uij(t)), 

u{t)c{t)u*{t) = diag( 7 i(t),7 2 (t), • • ■ ,7jv|-i(*)), (44) 
and introducing the new operators 

7Vf-l 

A M = E u?-(t)F if (45) 

on obtains the form ( )40l) of the generator. 

The structure ( HU1) takes into account the Hermiticity and trace preservation of the dy- 
namics, but does not guarantee its complete positivity. The formulation of necessary and 
sufficient conditions for the complete positivity of the dynamics of this generator is an im- 
portant unsolved problem. However, in the case that the rates are positive for all times, 

li(t) > 0, (46) 

the resulting dynamics is indeed completely positive, since the generator is then in Lindblad 
form for each fixed t > 0. 
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As an example let us consider the dynamical map given by Eqs. (l2"Tl) - (j2"4"l) . In this case 
the time-local generator takes the form [ijj] 



1 

£(t)p s = --S(t)[a+a-,p s ] 



+7(t) a-p s cr + 



-{a + a-,p s } , 



(47) 



where we have introduced the definitions 



7(t) = -23? 



G(t) 



G(t) 



) 



S(t) = -29 



G(t) 



G(t) 



) 



(48) 



With this generator Eq. (1371) represents an exact master equation of the model. The quantity 
S(t) plays the role of a time-dependent frequency shift, and j(t) can be interpreted as a 
time-dependent decay rate. Due to the time dependence of these quantities the process does 
generally not represent a dynamical semigroup, of course. Note that the generator is finite 
as long as G(t) 7^ 0; at the zeros of G(t) the inverse of the dynamical map $(i, 0) does not 
exist. An example is provided by the zeros of the function of Eq. (1291) which appear in the 
strong coupling regime 70 > A/2. 

We can also see explicitly how the standard Markov limit arises in this model. Considering 
the particular case (129|) . we observe that in the limit of small a = 70/A we may approximate 
G(t) ~ e -70 */ 2 . This approximation can also obtained directly from Eq. fl25l) by replacing 
the two-point correlation function f(t — ti) with the delta-function 7o5(t — ti), which is 
conventionally regarded as the Markovian limit. Equation ( T48l) then yields S(t) = and 
j(t) = 70, i.e., the generator (T47J) assumes the form of a Lindblad generator of a quantum 
dynamical semigroup. The quantity a can also be written as the ratio of the environmental 
correlations time te = A -1 and the relaxation time t r = j^ 1 of the system: 



Thus we see that the standard Markov condition (j34j) indeed leads to a Markovian semigroup 



a = — 



(49) 
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E. Divisibility of dynamical maps 



A family of dynamical maps $(£, 0) is defined to be divisible if for all i 2 > t\ > there 
exists a CPT map $(£2,^1) such that the relation 

$(t 2 ,0) = $(t 2 ,ti)$(ti,0). (50) 

holds. Note that in this equation the left-hand side as well as the second factor on the 
right-hand side are fixed by the given family of dynamical maps. Hence, Eq. (130]) requires 
the existence of a certain linear transformation $(£2, t\) which maps the states at time t\ to 
the states at time i 2 and represents a CPT map (which may be a unitary transformation). 
This definition differs slightly from the usual definition for the divisibility of a CPT map 
A, according to which A is said to be divisible if there exist CPT maps Ai and A 2 such 
that A = AiA 2 , where one requires that neither Ai nor A 2 is a unitary transformation, 
for otherwise the statement is trivial [33I] . There are many quantum processes which are 
not divisible. For instance, if $(ti,0) is not invertible, a linear map <&(i 2 ,ii) which fulfills 
Eq. ([50]) may not exist. Moreover, even if a linear map $(£ 2 ,£i) satisfying Eq. ( 150]) does 
exist, this map needs not be completely positive, and not even positive. 

The simplest example of a divisible quantum process is given by a dynamical semigroup. 
In fact, for a semigroup we have $(i, 0) = exp[£t] and, hence, Eq. (j50]) is trivially satisfied 
with the CPT map $(t 2 ,*i) = exp[£(t 2 - h)]. 

Consider now a quantum process given by the time-local master equation ([57]) with a 
time dependent generator (140]) . The dynamical maps can then be represented in terms of a 
time-ordered exponential, 



$(t,0) = Texp 



dt'K,{t') 







t>0, (51) 



where T denotes the chronological time-ordering operator. We can also define the maps 

-t 2 



Q(t2,h) = Texp 



dt'K(t') 



ti 



t 2 >h> 0, (52) 



such that the composition law $(£ 2 ,0) = $(^2; ti)&{ti, 0) holds by construction. The maps 
$(£ 2 ,ii) are completely positive, as is required by the divisibility condition (I50"]) . if and only 
if the rates 7i(i) of the generator (HO"]) are positive functions. Thus we see that divisibility is 
equivalent to positive rates in the time-local master equation Q. To prove this statement 
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suppose first that ji(t) > 0. The generator JC(t) is then in the Lindblad form for each fixed 
t > and, therefore, Eq. f|52|) represents completely positive maps. Conversely, assume that 
the maps defined by Eq. ([52]) are completely positive. It follows from this equation that the 
generator is given by 

d 



K{t) 



dr 



(53) 



T=0 



Since $(t + r, t) is completely positive for all t, r > and satisfies t) = I, this generator 
must be in Lindblad form for each fixed t, i.e., it must have the form fHOj) with 7i(t) > 0. 

The dynamical map $(t,0) given by Eqs. fl2Tl) - fl24l) is completely positive if and only if 
\G(t)\ < 1. It is easy to verify that 0) can be decomposed as in Eq. ([511 . where the map 
$(t 2 ,£i) is given by [3] 



Pooih) 

Pioih) 
poi(t 2 ) 



G(h) 



Poo(h 

G(t 2 ) 
G^) 
G*{t 2 



1 - 



G*{U 



Pio(ii), 
rPoi(*i)- 



(54) 



It follows from these equations that a necessary and sufficient condition for the complete 
positivity of $(t 2 ,^i) is given by 

\G{h)\ < |G(ti)|. (55) 

Thus we see that the dynamical map of the model is divisible if and only if |G(t)| is a 
monotonically decreasing function of time. Note that this statement holds true also for the 
case that G(t) vanishes at some finite time. The rate 7(t) in Eq. (j4~8]) can be written as 



7(*) 



d 



\G(t)\dt 



\G{t)\. 



(56) 



This shows that any increase of |G(t)| leads to a negative decay rate in the correspond- 
ing generator (j47p . and illustrates the equivalence of the non-divisibility of the dynamical 
map and the occurrence of a temporarily negative rate in the time-local master equation 
demonstrated above. 
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III. INFORMATION FLOW AND NON-MARKOVIAN QUANTUM DYNAMICS 



In the classical theory of stochastic processes a Markov process is defined by the Markov 
condition, which is a condition for the hierarchy of the n-point probability distribution func- 
tions pertaining to the process. Since such a hierarchy does not exist in quantum mechanics 
the question arises, how do quantum memory effects manifest themselves in the dynamical 
behavior of complex open quantum systems and how one can rigorously define and quantify 
non-Markovianity in the quantum case? Of course, such a definition must be independent of 



the specific mathematical representation of the open system's dynamics, i.e., it must 



De for- 



two 



mulated completely by means of the quantum dynamical map of the process. In Ref. 
different concepts of non-Markovianity have been proposed. The first one uses the fact that 
local quantum operations can never increase the entanglement between the open system and 
an isomorphic auxiliary system. By employing an appropriate measure for entanglement, a 
given dynamical evolution is then defined to be non-Markovian if a temporary increase of the 
entanglement measure takes place, and the size of this increase provides a measure for the 
degree of non-Markovianity. Within the second concept a quantum process is defined to be 
non-Markovian if and only if the dynamical map 0) is non-divisible. The corresponding 
measure for non-Markovianity quantifies the amount to which the dynamical map violates 
the divisibility condition, thus representing a measure for the non-divisibility character of 
the process. 

Here, we discuss a recent approach which defines and quantifies the emergence of quantum 
memory effects entirely in terms of the exchange of information between the open system 
and its environment 5|. In order to quantify this exchange of information we will use the 
trace distance as a measure for the distance between quantum states. 

A. Trace distance and distinguishability of quantum states 

The trace norm of a trace class operator A is defined by 

= tr|A| =WIU. (57) 
If A is self-adjoint with eigenvalues a;, this formula reduces to 

i 
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Hence, the trace norm of a self-adjoint operator is equal to the sum of the moduli of its 
eigenvalues (counting their multiplicities). The trace norm leads to a natural and useful 
measure for the distance between two quantum states represented by positive operators p 1 
and p 2 with unit trace, which is known as the trace distance: 

D{p\p 2 ) = \\\p^-p 2 \\ = l -tAp 1 - P \ (59) 

The trace distance is well defined and finite for all pairs of quantum states and provides a 
metric on the space S(7i) of physical states. We list some of the most important mathe- 
matical properties of this metric which will be needed later on (most of the proofs may be 



found in Refs. 



10 



oil)- 

1. The trace distance between any pair of states satisfies 

< D(p\p 2 ) < 1. (60) 

Of course, we have D(p l ,p 2 ) = if and only if p 1 = p 2 , while the upper bound is 
reached, i.e., D(p l ,p 2 ) = 1 if and only if p 1 and p 2 are orthogonal, which means that 
the supports of p 1 and p 2 are orthogonal. (The support is defined to be the space 
spanned by the eigenstates with nonzero eigenvalue.) Moreover, D(p l , p 2 ) is obviously 
symmetric in the input arguments, and satisfies the triangular inequality: 

D(p\p 2 )<D(p\p 3 )+D(p 3 ,p 2 ). (61) 

2. If p 1 = and p 2 = \i/j 2 ){i/j 2 \ are pure states the following explicit formula for 
the trace distance can easily be derived, 

D{p\p 2 ) = Vl-K^ 2 )! 2 - (62) 

If the underlying Hilbert space is two-dimensional (qubit), spanned by basis states 
|1) and |0), the trace distance between two states with matrix elements pi- and p?- is 
found to be 

D{p\p 2 ) = ^+W, (63) 

where a = p\ x — p\ x is the difference of the populations, and b = p\ Q — p 2 is the 
difference of the coherences of the two states. 
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3. The trace distance is sub-additive with respect to tensor products of states which 
means that 

D(p x ®cx\p 2 ® a 2 ) < D(p\p 2 ) + D{a\ a 2 ). (64) 

In addition we have 

D(p l ®a,p 2 ®a) =D{p\p 2 ). (65) 

4. The trace distance is invariant under unitary transformations U, 

D{Up l U\ Up 2 U^) = D(p\p 2 ). (66) 

More generally, all trace preserving and completely positive maps, i.e., all trace pre- 
serving quantum operations A are contractions of the trace distance, 

D{kp\kp 2 )<D(p\p 2 ). (67) 

Note that the condition of the trace preservation is important here, and that this 



inequality also holds for the larger class of trace preserving positive maps [341 ] . 

5. The trace distance can be represented as the maximum of a certain functional: 

D(p\p 2 ) = maxtr{U(p 1 -p 2 )}. (68) 

The maximum is taken over all projection operators EL Alternatively, one can take 
the maximum over all positive operators A with A < I. Note that this formula is 
symmetric, i.e., we can also write D(p l ,p 2 ) = maxntr{IT(p 2 — p 1 )} where, however, 
the maximum is then assumed for a different projection EL 

The trace distance between two quantum states p l and p 2 has a direct physical interpre- 
tation which is based on the representation ( l68l) . Consider two parties, Alice and Bob. Alice 
prepares a quantum system in one of two states p 1 or p 2 with probability 1 each, and then 
sends the system to Bob. It is Bob's task to find out by a single measurement on the system 
whether the system state was p 1 or p 2 . It turns out that Bob cannot always distinguish 
the states with certainty, but there is an optimal strategy which allows him to achieve the 
maximal possible success probability given by 

P^ = \[l + D{p\p 2 )]. (69) 
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Thus we see that the trace distance represents the bias in favor of the correct state identi- 
fication which Bob can achieve through an optimal strategy. The trace distance D(p l ,p 2 ) 
can therefore be interpreted as a measure for the distinguishability of the quantum states 
p 1 and p 2 [a-llQ]. 

According to Eq. (I68p Bob's optimal strategy consists in measuring the projection n for 
which the maximum in this relation is assumed, and in associating the outcome n = 1 
with the state p 1 , and the outcome II = with the state p 2 . Under the condition that 
the system state was p 1 he then has correctly identified this state with probability trjllp 1 }, 
while under the condition that the system state was p 2 his answer is correct with probability 
tr{(7 — n)p 2 }. Since both possibilities occur with a probability of ~ we obtain the success 
probability 

P m ax = ^tr{np 1 } + itr{(/-n)p 2 } 

^[i + Mn^-p 2 )}] 

= 1 -[l + D(p\p 2 )], (70) 

which proves Eq. ( 1691) . We further see that a state identification with certainty, P max = 1, 
can be achieved if and only if D{p 1 ,p 2 ) = 1, i.e., if and only if p l and p 2 are orthogonal, in 
which case Bob's optimal strategy is to measure the projection II onto the support of p 1 (or 
ofp 2 ). 



B. Definition of non-Markovian quantum dynamics 

As we have seen in Sec. IIII Al the trace distance Dip 1 , p 2 ) between two quantum states p 1 
and p 2 can be interpreted as the distinguishability of these states. An important conclusion 
from this interpretation is that according to the contraction property of Eq. (|67|) no com- 
pletely positive and trace preserving quantum operation can increase the distinguishability 
between quantum states. Consider two initial states p^(0) and p|(0) of an open system and 
the corresponding time evolved states 

p 1 f(t) = $(t,0)tf(0). (71) 

Since the dynamical maps $(t, 0) are completely positive and trace preserving, the trace 
distance between the time evolved states can never be larger than the trace distance between 
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the initial states, 

D( P m />!(*)) <£>(^(0), 4(0)). (72) 

Thus, no quantum process describable by a family of CPT dynamical maps can ever increase 
the distinguishability of a pair of states over its initial value. Of course, this general feature 
does not imply that D(p 1 s (t), p 2 s (t)) is a monotonically decreasing function of time. 

We can interpret the dynamical change of the distinguishability of the states of an open 
system in terms of a flow of information between the system and its environment. When 
a quantum process reduces the distinguishability of states, information is flowing from the 
system to the environment. Correspondingly, an increase of the distinguishability signifies 
that information flows from the environment back to the system. The invariance under 
unitary transformations fl66l) indicates that information is preserved under the dynamics of 
closed systems. On the other hand, the contraction property of Eq. (|67|) guarantees that the 
maximal amount of information the system can recover from the environment is the amount 
of information earlier flowed out the system. 

Our definition for quantum non-Markovianity is based on the idea that for Markovian 
processes any two quantum states become less and less distinguishable under the dynamics, 
leading to a perpetual loss of information into the environment. Quantum memory effect 
thus arise if there is a temporal flow of information from the environment to the system. 
The information flowing back from the environment allows the earlier open system states to 
have an effect on the later dynamics of the system, which implies the emergence of memory 
effects 

In view of these considerations, we thus define a quantum process described in terms of 
a family of quantum dynamical maps $(t, 0) to be non-Markovian if and only if there is a 
pair of initial states Ps' 2 (0) such that the trace distance between the corresponding states 
Ps 2 {t) increases at a certain time t > 0: 

a(t,pf(0))^jD(p 1 s (t),p 2 s (t))>0, (73) 

where a(t, Ps' 2 (0)) denotes the rate of change of the trace distance at time t corresponding 
to the initial pair of states. 

This definition for quantum non-Markovianity has many important consequences for the 
general classification of quantum processes. In particular, it implies that all divisible families 
of dynamical maps are Markovian, including of course the class of quantum dynamical 
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semigroups. To prove this statement suppose that $(£, 0) is divisible. For any pair of initial 



states Pg 2 (0) we then have 

pf{t + r) = $(t + r, t)pf(t), t, r > 0. (74) 

Since <&(i + t, t) is a CPT map we can apply the contraction property (|67|) to obtain: 

D(pl(t + r),p 2 s (t + r)) <D(pl(t),pUt)). (75) 

This shows that for all divisible dynamical maps the trace distance decreases monotonically 
and that, therefore, these processes are Markovian. 

Hence, we see that non-Markovian quantum processes must necessarily be described 
by non-divisible dynamical maps and by time-local master equations whose generator ( 140 p 
involves at least one temporarily negative rate 7i(t). However, it is important to note that 
the converse of this statement is not true, i.e., there exist non-divisible dynamical processes 
which are Markovian in the sense of the above definition. For such processes the effect 
of the contribution of the decay channels with a negative rate is overcompensated by the 
channels with a positive decay rate, resulting in a total information flow directed from the 
open system to the environment. Further implications are discussed in 35j, and specific 



examples for this class of processes are constructed in 3y, l37l |. 

We have interpreted the change of the distinguishability of quantum states as arising 
from an exchange of information between the open system and its environment. If the 
distinguishability decreases, for example, we say that information is flowing from the open 
system into the environment. The more precise meaning of this statement is that the avail- 
able relative information on the considered pair of states, that is, the information which 
enables one to distinguish these states is lost from the open system, i.e., when only measure- 
ments on the open system's degrees of freedom can be performed. This does not imply that 
the information lost is now contained completely in the reduced states of the environment, 
but instead, this information can also be carried by the system-environment correlations. 
And vice versa, if we say that information flows from the environment back to the system, 
this means that the distinguishability of the open system states increases because informa- 
tion is regained by the open system, information which was carried by the environmental 
degrees of freedom or by the correlations between the degrees of freedom of system and 
environment (see also Sec. HVj) . 
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C. Construction of a measure for the degree of non-Markovianity 



For a non-Markovian process described by a family of dynamical maps <&(£, 0) information 
must flow from the environment to the system for some interval of time and thus we must 
have a > for this time interval. A measure of non-Markovianity should measure the total 
increase of the distinguishability over the whole time evolution, i.e., the total amount of 
information flowing from the environment back to the system. This suggests defining a 
measure J\f(<$>) for the non-Markovianity of a quantum process through [5] 

W($) = max / dta^p 1 /^)). (76) 

The time integration is extended over all time intervals (a*, hi) in which a is positive and 
the maximum is taken over all pairs of initial states. Due to Eq. (173)) the measure can be 
written as 

M{$) = max :J2 M4(M,pI(&0) " ^(pM^pIK))]- (77) 

p ' (0) — ' 

i 

To calculate this quantity one first determines for any pair of initial states the total growth 
of the trace distance over each time interval (oj, h) and sums up the contribution of all 
intervals. J\f(<&) is then obtained by determining the maximum over all pairs of initial 
states. 

To illustrate our definition of non-Markovianity and the measure ( ITS"]) we again refer to 
the example of the dynamical map given by Eqs. ff2T]) -f l24|) . According to Eq. (163]) the time 
evolution of the trace distance corresponding to any pair of initial states Pg(0) and p|(0) is 
given by 

D( P Ut),P 2 s(t)) = \G(t)\y/\G(t)\W + \b\*, (78) 

where a = Pn(0) — Pn(0) and b = p\ (0) — Pi (0). The time derivative of this expression 
yields 

°(t,Ps (0))= - 7 __ 5i |G(i)|. (79) 

We conclude from this equation that the trace distance increases at time t if and only if the 
function |G(t)| increases at this point of time. It follows that the process is non-Markovian, 
A/'( < I ) ) > 0, if and only if the dynamical map is non-divisible, which in turn is equivalent to 
a temporarily negative rate j(t) in the time-local master equation with generator (I47p . 
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Consider, for example, the case of an exponential correlation function which leads to the 
expression (129]) for the function G(t). For small couplings, 70 < A/2, this function decreases 
monotonically. The dynamical map is then divisible, the rate 7(t) is always positive, and the 
process is Markovian. However, in the strong coupling regime, 70 > A/2, the function 
starts to oscillate, showing a non-monotonic behavior. Consequently, the dynamical map is 
then no longer divisible and the process is non-Markovian. Note that in the strong coupling 
regime the rate 7(t) diverges at the zeros of G(t). However, the master equation can still 
be used to describe the evolution between successive zeros and, therefore, the connection 
between non-Markovianity, divisibility and a negative rate in the master equation remains 
valid. There is thus a threshold 70 = A/2 for the system-reservoir coupling below which 
A/"($) = 0. One finds, as expected, that for 70 > A/2 the non-Markoviantiy increases 
monotonically with increasing system-environment coupling [7]. Moreover, it is easy to 
show with the help of Eq. (ITS"]) that the maximum in Eq. ( !76|) is attained for a = and 
\b\ = 1 [38]. This means that the optimal pairs of initial states correspond to antipodal 
points on the equator of the Bloch sphere representing the qubit. 

The non-Markovianity measure (ITSl) has been employed recently for the theoretical char- 
acterization and quantification of memory effects in various physical systems. An application 



to the information flow in the energy transfer dynamics o 



developed in 39[, and to ultracold atomic gases in 40[. Further applications to memory 



photosynthetic complexes has been 



effects in the dynamics of qubits coupled to spin chains |4JJ and to complex quantum sys- 



tems with regular and chaotic dynamics 42| have been reported. A series of examples and 
details of the relation between the classical and the quantum concepts of non-Markovianity 
are discussed in [43]. 

The non-Markovianity measure ( 1761) represents a novel experimentally accessible quantity. 



Quite recently, two experiments determining this quantity have 



jeen carried out employing 



photons moving in optical fibers or birefringent materials |44j. |45| . In these experiments the 
open system is provided by the polarization degrees of freedom of the photons, while the 
environment is given by their frequency (mode) degrees of freedom. It has been demonstrated 
that a complete determination of the measure is experimentally possible, including the 
realization of the maximization over initial states in the definition (T76j) . It was shown 
further that a careful preparation of the initial environmental state allows to control the 
information flow between the system and its environment and to observe the transition from 
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the Markovian to the non-Markovian regime through quantum state tomography carried out 
on the open system. This means that an experimental quantification and control of memory 
effects in open quantum systems is indeed feasible, which could be useful in the development 
of quantum memory and communication devices. 



IV. THE ROLE OF INITIAL SYSTEM-ENVIRONMENT CORRELATIONS 

As we have seen in Sec. Ill Bl the construction of a completely positive quantum dynamical 
map <3>(i, 0) acting on the open system's state space is based on the assumption of an un- 
correlated total initial state. In this section we will derive general inequalities which express 
the role of initial system-environment correlations in the subsequent dynamical evolution 



and discuss suitable strategies for the local detection of such correlations [46 ]. 



A. General relations describing the effect of initial correlations 

While in the general case of correlated initial states a quantum dynamical map can only 
be defined for a restricted set of states 47|, we can of course always introduce a linear map 

A(t, 0) : S(U S ® U E ) — ► S{U S ) (80) 

on the total system's state space S(Hs ®7~Le) by means of 

p(0) ^ ps (t) = A(t, 0)p(0) = tr E {U(t)p{0)U\t)} . (81) 

This map takes any initial state p(Q) of the total system to the corresponding reduced open 
system state pg(t) at time t. Since unitary transformations and partial traces are CPT 
maps, the composite maps A(£, 0) are again CPT maps. Considering a pair of total initial 
state p 1,2 (0) and the corresponding open system states at time t, 

ftf(t)=A(t, 0)^(0), (82) 

we then have by use of the contraction property (IBTj) for CPT maps: 

£(P^),PlW)< ^(0)^(0)), (83) 

which means that the distinguishability of the open system states can never be larger than 
the distinguishability of the total initial states. Subtracting the initial trace distance of the 
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open system states we obtain: 

£(P^),pfW)-^(Ps(0),p|(0)) 

<I(p\a),p\0))=D(p\0),p\0))-D(pi(0),pl(0)). (84) 

The increase of the trace distance between the open system states is thus bounded from 
above by the quantity 7(p 1 (0), p 2 (0)) on the right-hand side of this equation. This quantity 
represents the distinguishability of the total initial states minus the distinguishability of the 
initial open system states. It can be interpreted as the information on the total system states 
which is outside the open system, i.e., which is inaccessible through local measurements on 
the open system. Thus we see that the increase of the distinguishability of the open system 
states and, hence, the flow of information to the open system is bounded by the information 
which is outside the open system at the initial time. 

We now show how the upper bound of Eq. (184"]) is connected to the correlations in the 
initial states. To this end, we use twice the triangular inequality ( I6T1) to get: 

D( P \0),p 2 (0))<D( P \0),pl(0)®p E (0)) + D(p 2 (0),pl(0)®p E (0)) 

+D( P 1 s (0)®p 1 E (0),p 2 s (0)®P E (0)). (85) 
With the help of the subadditivity (JM|) of the trace distance we find 

D(p l s (0) ® p E (0),p 2 s (0) ® Pl(0)) < £>(pkO),p|(0)) + D(p E (0),p E (0)). (86) 

Combining this with Eqs. ( |85i) and ( |84l) we finally obtain 

D( P W),Ps(t)) ~ D(pl(0),p 2 s m < D(pH0),p' s (0) ® p E (0)) + £(p 2 (0),p|(0) ® p|(0)) 

+D(p E (0),p 2 E (0)). (87) 

For any total system state p the trace distance D(p,p$ ® Pe) between p and the product 
of its marginals ps ® Pe represents a measure for the correlations in the state p, which 
quantifies how well p can be distinguished from the corresponding product state ps <S> Pe- 
Thus, Eq. ( 1H71) demonstrates that a dynamical increase of the trace distance of the open 
system states over the initial value, 

£>G4(t),/4(t)) -D(p l s (0), pl(0)) > 0, (88) 
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implies that the corresponding initial environmental states p^ 2 (0) are different or that at 
least one of the total initial states p 1,2 (0) is correlated. 

As an example we consider the special case of an initial pair of states given by a correlated 
state p 1 (0) and the uncorrelated product of its marginals p 2 (0) = pg(0) ® Pe(0). Hence, we 
have Ps{0) = p|(0) and p^O) = p|(0), and Eq. (IHTJ) reduces to 

£>(Ps(*),pI(*)) < 23(^(0), p^(O) ® PKO)). (89) 

Thus, the increase of the trace distance between pg(t) and p|(t) (which is zero initially) is 
bounded from above by the amount of correlations in the initial state p 1 (0). 

We further remark that in the case of two uncorrelated states with the same environmental 
state, i.e., p l,2 {0) = Ps' 2 (0)®pe(0), the right-hand side of the inequality (IHTj) vanishes, and we 
are led again to the contraction property ( 1721) for CPT dynamical maps on the reduced state 
space. Thus, Eq. (IHTj) represents a generalization of the contraction property of dynamical 
maps on the reduced state space. 



B. Witnessing system-environment correlations by local operations 

ie local detection of cor- 



46] . First, one generates 



Equation ( 1871) leads to an experimentally realizable scheme for t 
relations in an unknown total system's initial state p x (0) as follows 
a second reference state p 2 (0) by applying a local trace-preserving quantum operation S, 

P 2 (0) = (£® l)p l (0). (90) 

The operation £ acts only on the variables of the open system, and can be realized, for 
example, by the measurement of an observable of the open system, or by a unitary trans- 
formation induced, e.g., through an external control field. It is easy to check that p x (0) and 
p 2 (0) lead to the same reduced environmental state, i.e., we have p E (0) = p E (0) and, hence, 
Eq. ([87D yields 

£(pM^pIW)-^(pM0),P 2 5 (0)) < J D(p 1 (0),p^(0)®p^(0))+ J D(p 2 (0),p|(0)®p|(0)). (91) 

This inequality shows that any dynamical increase of the trace distance between the open 
system states over the initial value implies the presence of correlations in the initial state 
p 1 (0). In fact, if one finds that the left-hand side of the inequality is greater than zero, then 
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either p 1 (0) or p 2 (0) must be correlated. If /^(O) was uncorrelated, then also p 2 (0) must be 
uncorrelated since it is obtained from p 1 (0) through application of a local operation. Thus, 
any increase of the trace distance of the open system states over the initial value witnesses 
correlations in p 1 (0). 

We note that this strategy for the local detection of initial correlations requires only local 
control and measurements of the open quantum system. It neither demands knowledge 
about the structure of the environment or of the system-environment interaction, nor a full 
knowledge of the initial system-environment state p 1 (0). Moreover, there is no principal 
restriction on the operation £ used to generate the second state p 2 (0), which makes the 



scheme very flexib 
reported recently 



e in practice. In fact, experimental realizations of the scheme have been 
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49]. Further examples and applications to the study of correlations 



in thermal equilibrium states are discussed in Ref . 50] • Moreover, by taking £ to be a pure 
dephasing operation the scheme enables the local detection of nonclassical correlations, i.e., 
of total initial states with nonzero quantum discord 511 ] . This fact has been shown very 
recently in Ref. |52j where also a statistical approach to initial correlations on the basis of 
random matrix theory has been developed. 



V. CONCLUSIONS 

We have discussed a definition for the non-Markovianity of quantum processes in open 
systems and developed a corresponding measures for the size of quantum memory effects. 
Our considerations are based on the quantification of the information flow between the open 
system and its environment in terms of the trace distance between quantum states of the 
open system. The great advantage of this distance measure is the fact that it admits a 
natural and clear physical interpretation as the distinguishability of the states through local 
measurements carried out on the open system. According to our definition the key feature 
of quantum non-Markovianity is the temporal increase of the distinguishability which can 
be interpreted as a backflow of information from the environment to the open system. As we 
have seen this concept allows a natural extension to the case of initial system-environment 
correlations. While a quantum dynamical map acting on the open system's state space does 
in general not exist in this case, the trace distance between pairs of states of the open system 
leads to a dynamical witness for the presence of initial correlations in the total system state. 
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Given that the initial correlations were created in the past from a product state through 
a system-environment interaction, the increase of the trace distance of the open system 
states over its initial value signifies that the open system regains information which was lost 
previously. 

There are of course several alternatives and possible modifications of the quantum mea- 
sure for non-Markovianity studied here. One possibility is to employ alternative distance 
measures for quantum states under which trace preserving quantum operations are con- 
tractive, such as the relative entropy or the Bures distance which is based on the fidelity 



10| | . In particular, the relative entropy represents a possibility which is natural both from 
an information theoretic perspective and from the point of view of nonequilibrium thermo- 
dynamics, since the negative rate of change of the entropy relative to an invariant thermal 



equilibrium state can be interpreted as entropy production 53|,|54J. However, a disadvantage 



of the relative entropy is given by the fact that it is often infinite, leading to singularities 
of the measure 5J. In several cases, in particular for infinite dimensional Hilbert spaces, 
the determination of the trace distance could be an extremely difficult task. An analytical 
formula for the trace distance is not even known for Gaussian quantum states. It seems 
that in those cases it is much easier to work with the Bures distance or the fidelity which 



also leads to useful lower and upper bounds for the trace distance [55(. The Hilbert-Schmidt 
distance which is technically much easier to deal with cannot be used for the quantifica- 
tion of non-Markovianity because trace preserving quantum operations are, in general, not 
contractive for this metric [56]. A further possibility is to define the measure for quantum 
non-Markovianity by means of alternative functionals of the family of dynamical maps which 
quantify the dynamical increase of the chosen distance measure for quantum states. 

Both the quantum measure for non-Markovianity and the witness for initial system- 
environment correlations studied here have been demonstrated very recently to be exper- 



imentally measurable quantities 
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45 



48 



49]. These experiments have paved the way 



for a series of further investigations of quantum memory effects in composite, multipartite 
open systems. A central goal of the quantum theory of open system is thus the design 
of appropriate schemes and models for their theoretical treatment and, in particular, the 
development of efficient numerical methods for the determination of quantum measures for 
non- Markovianity. 
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